Magnetic properties of Ag 2 VOP207: an unexpected spin dimer system 
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Magnetic properties of the silver vanadium phosphate Ag2VOP2C>7 are studied by means of 
magnetic susceptibility measurements and electronic structure calculations. In spite of the layered 
crystal structure suggesting ID or 2D magnetic behavior, this compound can be understood as a 
spin dimer system. The fit of the magnetic susceptibility indicates an intradimer interaction of 
about 30 K in perfect agreement with the computational results. Our study emphasizes the possible 
pitfalls in interpreting experimental data on structural basis only and points out the importance of 
microscopic models for the understanding of the magnetic properties of vanadium phosphates. 
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I. INTRODUCTION 

Low-dimensional spin-1/2 systems have been exten- 
sively studied during the last decade due to strong quan- 
tum fluctuations that may lead to unusual ground states 
and low-temperature properties. Magnetic frustration 
in low-dimensional spin systems additionally suppresses 
long-range ordering and results in fascinating phenomena 
such as the formation of a resonating valence bond (RVB) 
spin-liquid ground stated or the magnetoelectrical ef- 
fect in spin spirals^ Recent studies^ 6 -^ indicated the 
importance of combining experimental (measurement of 
thermodynamic properties, neutron and magnetic reso- 
nance studies) and computational (microscopic modeling 
via electronic structure calculations) approaches to gain 
an understanding of the unusual physical properties of 
low-dimensional spin-1/2 systems. 

Vanadium phosphates present a promising and mainly 
unexplored field for the search of novel low-dimensional 
(and, possibly, frustrated) spin systems. Structural re- 
ports on vanadium phosphates are often accompanied 
by magnetic susceptibility data- but few of the sys- 
tems gained thorough physical investigation, the no- 
table examples being (VO) 2 P 2 7 (Refs. [Hn© and 
VO(HP0 4 ) -0.5H 2 O3iii The careful study of these com- 
pounds revealed complex superexchange pathways and 
showed the failure of the straightforward interpretation 
of experimental results on the structural basis only i 10 i 14 

Recently, we studied magnetic properties of a novel 
vanadium phosphate Pb2VO(P04)2 and found an in- 
teresting realization of the frustrated square lattice sys- 
tem with both ferro- and antiferromagnetic interactions 
resulting in columnar antiferromagnetic ordering ! 15 ' 16 ! 17 
Thus, vanadium phosphates can reveal very unusual spin 
systems, and a detailed study of the respective com- 
pounds is of high interest. Below, we present the inves- 
tigation of the magnetic properties of a silver vanadium 
phosphate Ag 2 VOP 2 07.— We use both experimental and 
computational techniques in order to achieve a reliable 
description of the exchange couplings in the system un- 
der investigation. The results are discussed with respect 
to the structural correlations for magnetic interactions in 



vanadium phosphates. 

The outline of the paper is as follows. In Sec. UH we 
briefly describe the crystal structure and analyze the pos- 
sible pathways of superexchange interactions. Method- 
ological aspects are given in Sec. Mil Sec. IIVI deals with 
magnetic susceptibility data, while Sec. [V] presents the 
results of band structure calculations and estimates for 
the exchange integrals. Experimental and computational 
results are compared and discussed in Sec. I VII followed 
by our conclusions. 



II. CRYSTAL STRUCTURE 

The crystal structure of Ag 2 VOP 2 07 is formed by 
[VOP 2 07] layers separated by silver cations. The struc- 
ture has monoclinic symmetry (P2i/c, a = 7.739 A, 
b = 13.611 A, c = 6.294 A, @ = 99.0°, Z = 4). One unit 
cell contains two layers and, in particular, four vanadium 
atoms - two from each of the layers (Fig. [I]). Every 
layer includes VGe octahedra joined by PO4 tetrahedra. 
Vanadium has the oxidation state of +4 corresponding to 
the electronic configuration d 1 , i.e., S = 1/2*^ Magnetic 
measurements down to 77 K indicate paramagnetic be- 
havior of Ag 2 VOP 2 07,— but, as we will show in the fol- 
lowing, considerable deviation from the Curie law starts 
immediately below this temperature. 

To understand the possible magnetic interactions in 
Ag 2 VOP 2 07, one has to consider the local environment 
of vanadium atoms as well as the connections between 
vanadium polyhedra. The VOg octahedra are strongly 
distorted due to the formation of a short V-0 bond typi- 
cal for V+ 4 (seeRef. i). The distortion of the octahedron 
results in a non-degenerate orbital ground state with the 
half-filled d xy orbital lying in the equatorial plane of the 
octahedron,— i.e., almost parallel to the [VOP2O7] lay- 
ers. 

The orientation of the half-filled orbital implies that 
interlayer interactions should be negligibly small while 
the magnitude of in-layer interactions depends on the 
connections between the vanadium octahedra. Accord- 
ing to Fig. [H V-O-V superexchange is impossible in 



2 




FIG. 1: (Color online) Crystal structure of Ag 2 VOP 2 7 (up- 
per panels) and a schematic image of an alternating zigzag 
chain (bottom panel). The upper left panel shows a single 
[VOP2O7] layer, while the upper right panel presents the 
stacking of the layers. In the upper left panel, faces of the 
two V06 octahedra are removed in order to illustrate the ori- 
entation of short V-O bonds (shown as thicker lines). Solid, 
dashed, and dotted lines in the bottom panel denote J[, J", 
and J2, respectively. 

Ag2VOP207; therefore, all the superexchange pathways 
must include PO4 tetrahedra. Assuming that the in- 
teractions via two consecutive tetrahedra (i.e., the py- 
rophosphate P2O7 group) are weak, we derive three pos- 
sible couplings: J{, J" , and J2 indicating a ID spin sys- 
tem. One can consider such a system as an alternat- 
ing zigzag chain (see bottom panel of Fig. [T|) , since two 
nearest- neighbor interactions </{, J" are not equivalent 
by symmetry, while J2 is the next-nearest neighbor in- 
teraction. Interchain couplings are weak as they include 
P2O7 groups or single tetrahedra that do not match the 
half-filled orbital of vanadium. Below (Sec. [Vj, we justify 
the weakness of the interchain couplings by analyzing the 
local density approximation (LDA) band structure with a 
tight-binding model that accounts for all the interactions 
between nearest and next-nearest neighbors. 

III. METHODS 

Polycrystalline samples of Ag2VOP2C>7 were prepared 
according to the procedure reported in Ref. 18s. Stoi- 
chiometric amounts of AgNC>3, (NH4J2HPO4, and VO2 
were carefully mixed, placed in corundum crucibles and 
annealed in air. The regime of the annealing was as fol- 
lows: 12 hours at 200 °C, 1 day at 400 °C, and 1 day at 
500 °C with regrindings after each step. The final prod- 
uct was a single phase as checked by X-ray diffraction 
(XRD) (Huber G670f camera, CuK Q i radiation, Image- 
Plate detector). However, the XRD pattern revealed a 
rather high background signal that probably indicated 



low crystallinity of the prepared samples or the presence 
of amorphous impurities. Unfortunately, we failed to im- 
prove the quality of the samples, since further annealings 
(in air or in vacuum) did not lead to any notable changes 
in the XRD patterns or resulted in the appearance of non- 
identified impurity peaks. We also tried to prepare an 
isostructural compound Na2VOP2C>7 (Ref. [201 ) but the 
samples always contained a mixture of two polymorphic 
modifications [one is isostructural to Ag2VOP2C>7 while 
the other one has fresnoite-type structure (see Ref. |21|)]. 

Magnetic susceptibility was measured using Quantum 
Design SQUID between 2 and 300 K in the fields ^qH of 
0.5, 1, and 5 T. 

Scalar-relativistic density-functional (DFT) band 
structure calculations were performed using a full- 
potential local-orbital scheme^ (FPL05. 00-19) and 
the exchange-correlation potential of Perdew and 
Wang^ Ag(4s,4p,4d,5s,5p), V(3s, 3p, M, 4s, 4p), 
P(2s, 2p, 3s, 3p, 3d), and 0(2s,2p, 3d) orbitals were em- 
ployed as the basis set, while all lower-lying orbitals were 
treated as core states. The k mesh included 864 points 
within the first Brillouin zone (296 in its irreducible 
part). The convergence with respect to the number of k 
points was carefully checked. 

First, a LDA calculation was performed employing the 
crystallographic unit cell and the full symmetry of the 
crystal structure {P2\/c). Orbital states relevant for 
magnetic interactions were chosen on the basis of the 
LDA band structure and the respective bands were ana- 
lyzed using a tight-binding model. Such analysis pro- 
vided an estimate for the magnitudes of all nearest- 
neighbor and next-nearest-neighbor exchange couplings 
in Ag2VOP207. We found that the values of interest 
were J[, J", J2 (see Fig. [IJ, in agreement with the em- 
pirical structural considerations listed in Sec. [ill 

Next we turned to the LSDA+l/ technique in order to 
calculate total energies for several patterns of spin order- 
ing and to give an independent estimate of the leading 
exchange interactions in Ag2VOP207. Such an estimate 
requires doubling of the lattice translation in the c di- 
rection (see Section [Vj) . The doubling of the crystallo- 
graphic unit cell results in Z — 8 and 104 atoms making 
the calculations very time-consuming and less accurate. 
Therefore, we slightly simplified the crystal structure in 
our LSDA+J7 calculations. The simplification was justi- 
fied by the similarity of LDA band structures and tight- 
binding models for the crystallographic unit cell and the 
supercell (i.e., the modified cell). 

The neighboring layers in Ag2VOP207 are not trans- 
lationally equivalent, but they are equivalent by sym- 
metry, hence they reveal identical in-layer interactions. 
Interlayer interactions are very weak (see Sec. [Vj , there- 
fore the unit cell with one layer is actually sufficient for 
the present computational purpose. We use a PI super- 
cell with a s = a = 7.739 A, b s = 8.5069 A (~ 0.66), 
c s = 2c = 12.588 A, and (3 S = (3 = 99°. The supercell in- 
cludes one [VOP2O7] layer that exactly matches the layer 
in the real structure. However, the local environment of 
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silver atoms is changed, and the b s parameter is slightly 
elongated as compared to 6/2 in order to avoid unrea- 
sonably short Ag-0 distances. The modification of the 
crystal structure is a strong but well justified approxima- 
tion, since silver atoms do not take part in superexchangc 
pathways, while the geometry of the V-P-0 framework 
(crucial for magnetic interactions) is unchanged. These 
conclusions are supported by similar tight-binding fits of 
the relevant bands in LDA band structures for the crys- 
tallographic unit cell and the supercell. The variation of 
the respective hopping parameters does not exceed 5%, 
hence the error for J's is below 10% (i.e., within the typ- 
ical accuracy of LSDA+C7). 

Another problematic point of LSDA+?/ calculations 
deals with the choice of the Coulomb repulsion param- 
eter Ud [we use the notation Ud in order to distinguish 
this parameter of the computational method from the 
model parameter J7 e ff for the Coulomb repulsion in the 
effective one-band model (see Sec. IVj)]. Below, we use 
several physically reasonable values of Ud (namely, 4, 5, 
and 6 eV) in order to check the influence of Ud on the 
resulting exchange integrals. The exchange parameter 
of LSDA+t/ is fixed at J = 1 eV, as usually it has no 
significant influence on the results. 
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FIG. 2: (Color online) Magnetic susceptibility of Ag 2 VOP2C>7 
measured at fields fj,oH of 0.5 T (empty circles) and 5 T (filled 
circles). Solid and dashed lines show the fits with the model 
of S = 1/2 alternating chain according to the first and second 
lines of Table HI respectively. The dash-dotted line is the fit 
with the HTSE (Ref. [H) above 10 K. 



IV. EXPERIMENTAL RESULTS 

Magnetic susceptibility curves for Ag 2 VOP 2 Oy are 
shown in Fig. [2] The curves reveal a maximum at 
T m ax ~ 20 K and a fast decrease from T max down to 
5 K indicating spin-gap-like behavior. Below 5 K the 
susceptibility increases (at low fields) or remains nearly 
temperature-independent (at 5 T). The samples under 
investigation are powders, therefore the low-temperature 
part of the data is dominated by the contribution of 
paramagnetic impurities and defects saturated at high 
field. The susceptibility curves above 2.2 K do not show 
any anomalies that could denote long-range magnetic 
ordering^ A weak bend is present on the low-field curve 
below 2.2 K, but the present data (collected above 2 K) 
are not sufficient to reveal the origin of this bend. 

Thus, the susceptibility data for Ag 2 VOP 2 07 reveal 
low-dimensional spin-gap-like behavior consistent with 
the scheme of magnetic interactions proposed in Sec. [TJ 
The alternation of nearest-neighbor interactions J[ and 
J" causes the formation of the spin gap. The high- 
temperature part of the susceptibility (above 60 K) shows 
the Curie- Weiss type behavior, although the 1/x vs. T 
curves are nonlinear due to considerable temperature- 
independent contribution. We fit the data above 100 K 
as x = Xo + C/(T - 6>) and find X o = -7.9(1) • 10~ 4 
cmu/mol, C = 0.350(1) emu-K/mol, and 9 = 11.0(4) K. 
The Curie constant C corresponds to the effective mo- 
ment of 1.62 fiB that is slightly lower than the expected 
spin-only value of 1.73 /is- The discrepancy may be 
caused by the presence of nonmagnetic impurities in the 
samples under investigation, and this explanation is ad- 



ditionally supported by high diamagnetic contribution 
Xo as well as by the slightly decreased g value (see be- 
low). We suggest that the high diamagnetic contribution 
is caused by the preparation procedure, since the samples 
are heated in air, i.e., in oxidative conditions favoring the 
formation of diamagnetic V +5 . Unfortunately, annealing 
in vacuum did not result in single-phase samples. 

Now, we fit the experimental data with model expres- 
sions. Basically, one should use the model of the alter- 
nating zigzag chain, as suggested by structural consider- 
ations in Sec. |TTJ However, only high-temperature series 
expansion (HTSE) results are available for the suscepti- 
bility of this models We find a perfect fit of the exper- 
imental data above 10 K (see Fig. (2) with g = 1.92(1), 
J[ = 31.3(1) K, J'{ = 0.8(3) K, J 2 = 5.0(6) K, and 
Xo = —8.7(1) • 10 emu/mol (temperature-independent 
contribution). Note that we have assumed J" < J[: at 
the stage of fitting, this choice is, of course, arbitrary, 
and one may invert the assumption without any change 
in the results. One should not be surprised by the use of 
the HTSE for temperatures well below J[ , as the expan- 
sion under consideration is valid down to T w 0.25 J (see 
Ref. [HI). Moreover, the data below the susceptibility 
maximum (i.e., below 20 K) are crucial for the stability 
of the fitting and the obtaining of reasonable exchange 
couplings. 

Unfortunately, the HTSE is inapplicable in the low- 
temperature region. To fit the decrease of the suscepti- 
bility between 20 and 5 K, we have to turn to a simplified 
model. If one neglects J 2 , Ag 2 VOP 2 07 becomes an alter- 
nating chain system. The experimental data above 5 K 
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TABLE I: The results of the fitting of magnetic susceptibil- 
ity data with an alternating chain model [eq. JTJ]. xo is a 
temperature-independent contribution, d denotes the Curie- 
like Ci/T contribution, g is the (/-factor, Ji labels the leading 
in-chain interaction, and a is the alternation parameter. Fit 
1 and fit 2 are marked according to Fig. [2] 

Xo (emu/mol) d (emu-K/mol) g Ji (K) a 
Fit 1 -8.3(1) • 10" 4 0.0100(2) 1.86(1) 33.1(1) 0.12(2) 
Fit 2 -9.0(1) ■ 10" 4 1.92(1) 31.6(1) 0.31(1) 



are fitted with the expression: 

X = XO + Q./T + Xalt.chain, (1) 

where xo is the temperature-independent contribution, 
Ci/T accounts for the Curie- like paramagnetism of im- 
purities and defects, while Xait. chain is the susceptibility 
of an S = 1/2 alternating chain as given in Ref. [2a 

The term Xalt.chain includes three variable parameters: 
g, Ji, and the alternation coefficient a. Again, we as- 
sume J" < J[, then J i — J[ and a = J" / ' J[. We find 
a perfect fit of the experimental data above 5 K (fit 1 in 
Fig. (2) with the parameters listed in the first row of Ta- 
ble HI Below 5 K, the experimental data deviate from the 
model due to the deviation of the paramagnetic impurity 
contribution from the Curie law. Note that one may skip 
the Ci/T term and find a reasonable fit above 8 K (fit 
2 at Fig. [21 the parameters are listed in the second row 
of Table fl} , but the data between 5 and 8 K are poorly 
described. However, the second fit presents a somewhat 
better g value that is close to the expected range of 1.93- 
1.96 (see Refs. [27] and l28l). On the other hand, g may 
be decreased due to the presence of nonmagnetic impu- 
rities in the sample under investigation (as shown by the 
Curie- Weiss fit above) or due to the neglect of J2. 

Thus, the present experimental data leave some un- 
certainty in the a value. In our opinion, one can not 
neglect the Curie-like paramagnetic contributions while 
dealing with the low-temperature susceptibility of poly- 
crystalline samples. Therefore, we suggest that the first 
fit with the extremely low a is more realistic as confirmed 
by the HTSE fit [and band structure calculations (see the 
next section)]. Note that both fits reveal a considerable 
difference between J[ and J" (at least by a factor of 3). 
However, the scenario of J" <g; J[ is rather unexpected, 
since superexchange pathways of the two interactions are 
very similar. 



V. BAND STRUCTURE 

The plot of LDA density of states is shown in Fig. [3] 
At large, the band structure of Ag 2 VOP207 is simi- 
lar to that of other vanadium compounds (see, for in- 
stance, Refs. [l^ and [29|): the states below —3.5 eV 
are dominated by oxygen orbitals, while the states near 
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FIG. 3: (Color online) Total and atomic resolved LDA density 
of states for Ag2VOP20y. The Fermi level is at zero energy. 
The inset shows the blowup of the image near the Fermi level. 
The solid filling in the primary figure marks the contribution 
of all vanadium orbitals, while the solid filling in the inset 
corresponds to V 3d xy orbitals only. 

the Fermi level have mainly vanadium 3d character with 
minor oxygen contribution. However, a set of narrow 
bands between —1 and —3 eV is a distinctive feature 
of Ag2VOP2C>7. These bands correspond to 4d orbitals 
of silver and are completely filled as one would readily 
expect for Ag +1 (4d 10 ). The LDA energy spectrum re- 
veals a narrow gap of about 0.015 eV that is definitely 
too small to account for the green color of Ag 2 VOP207. 
The strong underestimate (or even the lack) of the en- 
ergy gap is a typical failure of LDA in transition metal 
compounds due to an improper description of the cor- 
relation effects in the 3d shell (especially for half-filled 
orbitals, such as 3d xy in V +4 ). Realistic insulating spec- 
tra with E g — 2.0 — 2.3 eV are readily obtained by means 
of LSDA+C/ (see below). 

Below, we analyze the band structure in order to find 
the values of exchange couplings and to construct a re- 
liable model of magnetic interactions in Ag2VOP207. 
First, we focus on the LDA bands near the Fermi level 
and extract transfer integrals for vanadium orbitals using 
a tight-binding model. The transfer integrals (t) are in- 
troduced to the extended Hubbard model, and the corre- 
lation effects are taken into account explicitly via an effec- 
tive on-site repulsion potential U c ff . In the strongly cor- 
related limit t <C [/ ff, the half-filled Hubbard model can 
be reduced to a Heisenberg model for the low-lying ex- 
citations. Thus, we estimate magnitudes of all exchange 
interactions between nearest and next-nearest neighbors. 
Next, we give an independent estimate of several ex- 
change coupling constants by introducing the correlation 
effects in a mean-field approximation to the band struc- 
ture (LSDA+C7) and calculating total energies for differ- 
ent patterns of spin ordering. 

The four bands near the Fermi level have V d xv char- 
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TABLE III: Exchange integrals calculated with LSDA+t/" 



FIG. 



4: (Color online) LDA band structure of Ag 2 VOP2C>7 
near the Fermi level. Thick (green) lines show the 
fit of the tight-binding model. The Fermi level is at 
zero energy. The notation of k points is as follows: 
F(0,0,0), B(0.5,0,0), A(0.5, 0.5,0), Y(0, 0.5,0), Z(0, 0,0.5), 
D(0.5, 0,0.5), E(0.5,0.5,0.5), C(0,0.5,0.5) (the coordinates 
are given along k x , k y , and k z axes in units of the respec- 
tive reciprocal lattice parameters). 



TABLE II: Leading hopping parameters of the tight-binding 



model and the resulting values of J 
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actcr (see the inset Fig. [3]), in agreement with the 
qualitative discussion in Sec. HH We fit a four-band 
tight-binding model to these bands (Fig. Ql and use 
the resulting hopping parameters to estimate antifer- 
romagnetic contributions to the exchange integrals as 
jAfm _ 4£2/jj eB , jj cS i s the effective on-site Coulomb 
repulsion. According to Ref. IS, U e s lies in the range of 
4 — 5 eV for vanadium oxides, and we use U c g = 4.5 eV in 
the following discussion as a representative value. Note 
that the change of U c g results in a simple scaling of the 
exchange integrals, therefore the ratios of J's remain con- 
stant. 

The tight-binding model included seven in-layer 
hoppings (i.e., all nearest-neighbor and next-nearest- 
neighbor hoppings) and an interlayer hopping t±. The 
resulting t± « 0.001 eV is very small, as one can readily 
see from Fig. SJ the bands are close to double degeneracy 
since two vanadium atoms of one layer weakly interact 
with the respective atoms of the neighboring layer. The 
leading in-layer hoppings t[, t", and t 2 (see Fig. Q} arc 
listed in Table HIl while other in-layer hoppings are small 
(below 0.008 eV) and can be neglected. Thus, the anal- 
ysis of the LDA band structure confirms the empirical 
structural considerations listed in Sec. [Ill 

The leading transfer integral t[ corresponds to 

j|AFM „ 37 K whUe j,/AFM and jAFM &re weaker by 

an order of magnitude. This result is in a remarkable 
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agreement with the experimental data matching both the 
value of the strongest coupling and the striking difference 
between nearest-neighbor interactions J[ and J". Note, 
however, that the tight-binding model does not allow us 
to distinguish between two structurally similar superex- 
change pathways, and one can also fit the bands with 
t" ^> t[. Now, we proceed to a different computational 
technique that helps us to resolve t[ , t'[ and provides ad- 
ditional confirmation for the difference of the respective 
exchange integrals. 

The hopping parameters (indicating the width of the 
half-filled bands) are two orders of magnitude smaller 
than the Coulomb repulsion; therefore, we may use the 
LSDA+J7 technique to account for correlation effects 
within band structure calculations in a mean-field ap- 
proximation. Four patterns of spin ordering (see Fig. [5]) 
were used in order to estimate Jj, J", and J 2 . The re- 
sulting exchange integrals are listed in Table IIIIl The 
absolute values of J's depend on Ud considerably, but 
qualitatively the results are similar. We find that J[ is 
indeed the strongest coupling, it exceeds J", Ji at least 
by a factor of 6. J{ is the interaction corresponding to 
longer V-V separation (5.29 A vs. 4.92 A for J") as 
shown in Fig. [T] The small value of J 2 is consistent with 
the tight-binding results and provides a further justifica- 
tion for the fitting of experimental data with the alter- 
nating chain model (see Sec. IIV[) . All the interactions in 
Ag2VOP2C>7 are antiferromagnetic. 

As we have mentioned in Sec. IIIIl the choice of Ud 
is a problematic point of the LSDA+J7 approach. The 
Ud parameter is usually fitted to experimentally observed 
properties (e.g., energy gaps, magnetic moments) or es- 
timated by means of constrained LDA calculation (see, 
for example, Ref. [29h . Basically, the choice of Ud re- 
mains somewhat ambiguous and depends on the partic- 
ular experimentally observed quantities used for the fit- 
ting as well as structural features and the computational 
method. 

The results for Ud = 6 eV show the best agreement 
with the experimental and tight-binding estimates of J's. 
Note also that the largest energy gap of 2.3 eV found for 
Ud = 6 eV seems to be consistent with the green color of 
Ag2VOP207. These findings are somewhat unexpected 
since we used U e g — 4.5 eV, and even smaller Ud — 3.6 eV 
had been applied in LSDA+[/ calculations for vanadium 
oxides before.-^ However, one should always keep in mind 
that C/ e fi and Ud are different quantities having different 
physical meaning. U e s is an effective repulsion in a mixed 
vanadium-oxygen band, while the repulsive potential Ud 
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FIG. 5: (Color online) Four patterns of spin ordering em- 
ployed for the calculation of J's with LSDA+17. Dash-dotted 
rectangles show the repetition unit (i.e. the supercell). Circles 
denote vanadium atoms; full and empty circles mark different 
spin directions. 



is applied to V 3d orbitals only. Thus, Ud should exceed 
U c ff, since U e s corresponds to "molecular" V-0 orbitals 
having larger spatial extension as compared to that of 
atomic V 3d orbitals. As for the previous LSDA+f/ cal- 
culations for vanadium oxides^ the difference in com- 
putational methods and crystal structures may be rel- 
evant. The linearized muffin-tin orbital (LMTO) basis 
is used in Ref. 29, while we apply a different basis of 
local orbitals. Moreover, Korotin et al^ studied lay- 
ered compounds containing V +4 in square-pyramidal co- 
ordination, while vanadium is surrounded by six oxygen 
atoms in Ag2VOP20y. Thus, the difference in the op- 
timal Ud value is plausible, although the origin of this 
difference is not completely clear. 



VI. DISCUSSION 

Experimental data and band structure calculations 
provide a consistent microscopic scenario for the mag- 
netic interactions in Ag2VOP2C>7. One of the nearest- 
neighbor couplings (J{) is about 30 K and exceeds 
the other couplings at least by a factor of 6. Thus, 
Ag2VOP207 may be considered as a system of weakly 
coupled dimers. Indeed, a simple model of isolated 
dimers [with additional terms xo + Ci/T similar to eq. 
(JTJ] with g — 1.83 and J = 33 K perfectly fits the ex- 
perimental data above 5 K. This fit completely matches 
fit 1 in Fig. O Interdimer interactions J" and J2 may 
be responsible for the slight variation of the g value, and 
they certainly influence the possible onset of the long- 
range ordering. The HTSE fit and band structure cal- 
culations indicate that both J" and J2 are antiferromag- 
netic, therefore frustration should be expected. How- 
ever, it will not change the main feature of Ag2VOP2 0y 
- a spin-gap-type magnetic behavior due to strong in- 
tradimer coupling. 

The identification of weakly coupled spin dimers in 
Ag 2 VOP207 is surprising, since two nearest-neighbor in- 
teractions - J{ and J" — have very similar superexchange 
pathways. The differences in interatomic distances and 



angles for these V-O-P-O-V connections do not exceed 

0. 05 A and 10°, respectively. Nevertheless, the tiny struc- 
tural changes strongly influence the values of the ex- 
change integrals. 

As we have mentioned in the introduction, superex- 
change pathways in vanadium phosphates are sometimes 
non-straightforward. However, some clarity may be 
achieved by considering two assumptions. First, V-O-P- 
O-V superexchange pathways are often much more effi- 
cient than V-O-V ones. Second, the magnitude of inter- 
actions involving PO4 tetrahedra can be estimated using 
structural parameters [via the so-called magnetostruc- 
tural correlations, see Refs. and [3l|. These corre- 
lations suggest that double phosphate bridges (i.e., two 
PO4 tetrahedra linking two VOg octahedra) are usually 
more efficient compared to single bridges, while the inter- 
actions via double bridges depend on the relative orienta- 
tion of vanadium and phosphorous polyhedra quantified 
by several geometrical parameters. The present study 
may provide a good test for such structural considera- 
tions concerning magnetic interactions in systems with 
complex superexchange pathways. 

In case of Ag 2 VOP207, double phosphate bridges are 
relevant for both J[ and J". Most geometrical parame- 
ters are similar for these interactions, and the only differ- 
ence deals with the in-plane displacement of the VOg oc- 
tahedra (plane implies the equatorial plane of octahedra, 

1. e., the plane of the half-filled d xy orbital). Such dis- 
placement favors J[ rather than J" , although one would 
not expect a considerable difference between the two in- 
tegrals. However, our experimental and computational 
results suggest the scenario of J[/ J" > 6 resulting in a 
system of weakly coupled dimers rather than in an al- 
ternating chain. Basically, the structural considerations 
suggest correct trends, but quantitative estimates require 
more sophisticated techniques. 

In this paper, we tried to use an adequate computa- 
tional approach to the study of vanadium phosphates. 
All the computational results available are based on very 
simplified techniques such as extended Hucke l 30 i 32 ' 33 or 
DFT GGA calculations 3 ^ 3 ^ for small clusters (typi- 
cally, spin dimers - two VOg octahedra and bridging 
PO4 tetrahedra). Now, we show that vanadium phos- 
phates are readily accessible to full-potential band struc- 
ture calculations that are known as a very efficient tool 
for the study of magnetic properties of transition metal 
oxides ^iLl i 19 ! 29 ! 36 ! 37 Full-potential calculations provide 
reasonable accuracy and set the computational results 
on much safer grounds compared to cluster calculations, 
since the crude cluster approximation is avoided and cor- 
relation effects (which are very important for the prop- 
erties of transition metal compounds) can be included. 
We are looking forward to the further application of full- 
potential calculations to vanadium phosphates and to the 
discovery of novel unusual spin systems in these com- 
pounds. 

It is worth noting that Ag2VOP2 07 may be an inter- 
esting realization of the alternating zigzag chain model. 
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This model, also known as dimerized zigzag chain or 
dimerized frustrated chain, has been extensively stud- 
ied theoretically with a particular interest to its excita- 
tion spectrum i 38 ' 39 i 40 ' 41 i 42 ' 43 However, the experimental 
access to the model is still scarce and basically limited 
by the spin Peierls phase of CuGeC>3 that corresponds 
to the weakly dimerized case. As we have stated above, 
Ag2VOP2C>7 is strongly dimerized, and the dominating 
spin-gap properties are observed. Nevertheless, the frus- 
tration of weak interdimer couplings may be manifested 
in high-field properties of the material or its magnetic 
excitations. 

Using the representative values of J[ ~ 30 K, J" ~ 
J2 ~ 5 K, we find a/ ~ 0.17, S ~ 0.7 and estimate 
the position of Ag2VOP2C>7 on the phase diagram of 
the model [a/ = Jij J\ is the frustration ratio, while 
S = (J' - J") /(J' + J") is the dimerization ratio]. Ac- 
cording to Refs. |4J andl45l. this position implies the pres- 
ence of magnetization plateau, while Ref. |40j suggests the 
formation of bound singlet and triplet states above the 
elementary triplet excitation that corresponds to the spin 
gap. Thus, a detailed study of the high-field properties 
of Ag2VOP2C>7 and neutron scattering experiments prob- 
ing the magnetic excitations of this compound could be 
interesting. 

In conclusion, we studied magnetic properties of 
Ag2VOP2C>7 using magnetic susceptibility measurements 



and full-potential band structure calculations. We reach 
a perfect consistency between experimental and compu- 
tational results and describe Ag 2 VOP207 as a system of 
weakly coupled dimers on the alternating zigzag chain, al- 
though the strong dimerization is rather surprising from 
the structural point of view. The frustration of inter- 
dimer interactions may lead to a number of unusual prop- 
erties that deserve a further experimental study. We 
present the first application of full-potential band struc- 
ture calculations to vanadium phosphates. The calcu- 
lations show reasonably accurate and highly reliable re- 
sults, hence they should be an efficient tool in the study 
of related compounds. 
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